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For a set of TV identical massive boson wavepackets with 
optimal initial quantum mechanical localization, we calcu- 
late the Hanbury-Brown/Twiss (HBT) two-particle correla- 
tion function. Our result provides an algorithm for calculat- 
. ing one-particle spectra and two-particle correlations from an 
Cllrbitrary phase space occupation (qj, Pi,t»)i=i,jv as e.g. re- 
turned by event generators. It is a microscopic derivation of 
One result of the coherent state formalism, providing explicit 
finite multiplicity corrections. Both the one- and two-particle 
jgbectra depend explicitly on the initial wavepacket width a 
^wnich parametrizes the quantum mechanical wavepacket lo- 
"^salization. They provide upper and lower bounds which sug- 
,_gest that a realistic value for a has the order of the Compton 
^"^avelength. 

<^f}ACS numbers: 25.75.+r, 07.60.1y, 52.60.+h 
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Two-particle correlations C(Q, K) of identical parti- 

ies are the only known observables giving access to the 

pace-time structure of the particle emitting source in 

leavy ion collisions. Their interpretation is based on the 

O^sult of the coherent state formalism jl]]^] which reads 

"^5a jthe plane wave approximation for a large number of 
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this setting, an Hanbury-Brown/Twiss (HBT) inter- 
^rometric analysis aims at extracting from the correla- 
tor C(Q,K) as much information as possible about the 
space-time emission function S(x,K). Since this emis- 
sion function cannot be reconstructed unambiguously 
from C(Q, K) ||, (1) is mainly used in the study of model 
emission functions S(x,K). These studies have clarified 
to a considerable extent the question which geometri- 
cal and dynamical source characteristics are reflected in 
which particular momentum dependencies of the corre- 
lator (cf. m and refs. therein). A comparison with mea- 
sured correlations then allows to constrain the class of 
source models consistent with data. 

Microscopic event generators are one important tool 
to generate model emission functions. Here, we do 
not discuss in how far existing event generators (e.g. 
RQMD §}, VENUS f|, ARC @) provide an internally 
consistent calculation of the phase space distribution. 
None of them propagates (anti)-symmetrized N-particle 



states from first principles, and the resulting difficulties 
in calculating 2-particle correlations have been discussed 
recently in great detail [Q. The typical event generator 
output is a set £ of phase space points at given times 
z i = (QijPiji-t) which one associates with the "points of 
last interactions" . However, the Heisenberg uncertainty 
principle allows to interpret the z, only as mean positions 
of boson wave packets. To specify the localization of 
these wavepackets in phase space, at least one additional 
parameter is needed, e.g. the initial spatial wavepacket 
width a. Irrespective of how the phase space occupation 
has been obtained, we shall take the set £ and the width 
a as initial condition for the present investigation: £ and 
a define the boson emitting source. For notational sim- 
plicity, we restrict our discussion to one particle species, 
negative pions, say. 

The problem in associating an emission function 
S(x, K) to the distribution £ is that £ is a discrete phase 
space distribution of on-shell particles. In contrast, the 
emission function S(x, K) of the coherent state formal- 
ism is a continuous distribution which allows for off-shell 
momenta K. Often, one circumvents this problem by 
approximating the Pratt algorithm j|] by an ad hoc pre- 
scription: each particle pair (i,j) is weighted with a prob- 
ability pij, qi being 4- vectors % = (tf,q<), 

C(AQ, AK) = N(A l AK} Pii , (2a) 

Pij = 1 + cos(( Pi - pj) ■ ( qi - qj)) . (2b) 

Here, C(AQ,AK) denotes the 2-particle correlator for 
pairs whose relative and average pair momenta pi — pj, 
|(Pi + Pj) lie in the bin AQ, AK. N(AQ, AK) is the 
corresponding number of particle pairs. A tentative ar- 
gument to justify the prescription (|J) is that p^j coin- 
cides with the formal Born probability density of 
the Bose-Einstein symmetrized 2-particle plane wave 

Pi] = ^*{Qi,Qj,Pi,Pj)^{li,Qj,Pi,Pj) , (3a) 
*(<ft, qj,Pi,Pj) = -L(e*«*+<W* + e<w«+*«") . (3b) 

However, the prescription (^) based on the ansatz (||) is 
inconsistent || with the result (1) of the coherent state 
formalism: The correlator in (1) is always larger than 
unity ||]. In contrast, the expression @ can drop below 
unity in the region of sufficiently large relative momenta 
||. Also, the prescription (||) is difficult to reconcile with 
quantum mechanical localization requirements since the 
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plane wave (3b) cannot be an eigenstate for both the 
position and momentum operator. 

In what follows, we take the quantum mechanical lo- 
calization of bosons into account by associating to the 
phase space emission points z; free Gaussian wavepack- 
ets of initial spatial width a, [10 11 
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This one boson state ( pa] ) is optimally localized around 
(qi, Pi) in the sense that it saturates the Heisenberg un- 
certainty relation Ax ■ Ap x = 1, with Axi = a at time 
t — ti. The time evolution of (Q) is the free unperturbed 
evolution determined by the Hamiltonian Hn — -s—, A 
being the Laplacian. Since the z-th and j-th boson are 
identical, we associate to the two emission points z$ and 
Zj the symmetrized two boson wave function $^ (X, Y, t) 
(the normalization factor is omitted and plays no role in 
what follows) 

* ii (X,Y,«)=/WpC,t)/W(Y,*) 

+/W(Y,t)/W(X,t). (5) 

We now derive an algorithm for calculating one-particle 
spectra v(P) and two-particle correlations C(Pi,P2) 
from an arbitrary initial phase space distribution £ of 
best localized boson wavepackets /if . Our first step is 
to calculate for two identical bosons the detection proba- 
bility at time t at the positions X and Y with momenta 
Pi, P2 respectively. This is given by the two-particle 
Wigner phase space density Q 

Wy(X,Y,Pi,P 2 ,t) ^^ 3 (X,Y,t)(27T) 6 S ( - 3 \P 1 - Pi) 
x<5( 3 )(P 2 -P 2 )$* 3 (X,Y,i) 

d 3 Xi d 3 Yi$(X 
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The corresponding probability to detect these bosons 
with momenta Pi and P 2 irrespective of their position 
is 

P ij (P 1 ,P 2 ) = J d 3 Xd 3 YW lJ (X,Y,P 1 ,P 2 ,t) 

= w i (P 1 ,P 1 )w j (P2,P 2 ) 
+iu i (P 2 ,P 2 )w J (P 1 ,P 1 ) 

+2ifl i (Pi,Pa)u> J -(Pi,P 2 ) 

x cos ((qi-qj)- (P% - P 2 )) , (7a) 

Wl (Pi,P 2 ) =e-^( p i - p 2) 2 Si ( K ), (7b) 
Sl (K) = 2 3 (w 2 )ie- CT2 ^- K ) 2 . (7c) 



Here, denotes 4-vectors Pi — (^PfjPj). We note 
that Pij is independent of the detection time t, i.e., only 
the correlations which exist already at emission are mea- 
sured at time t in the detector. This ^-independence is 
a consequence of the free time evolution; it is lost if fi- 
nal state interactions are included in the evolution of the 
wavepackets (^). Neglecting higher order symmetriza- 
tions, we define the (unnormalized) two pion correlation 
i?(Pi,P 2 ) for a set of N phase space points Z{ by sum- 
ming the probabilities P^ over all ^N(N — 1) pairs 

i?(Pl,P 2 ) =^Py(P 1 ,P 2 ). 

= i/(Pi)i/(P 2 )-2T c (Pi,P 2 ) 



N 



+ |^ U ; 4 (Pi,P 2 )e^( f!l - B2 )- iq '( Pl - p2) 

i=l 
N 

I/(P)=5>(P). 
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Here, Si(P) is the one-particle probability ( f7c| ) that a 
boson in the state /jf-* is detected with momentum P. 
Accordingly, v(P) is the one-particle spectrum of the dis- 
tribution £ with spatial localization a. The contribution 
T c to i?(Pi,P 2 ) corrects for the fact that the sums in 
the other two terms of (Q) include the N identical pairs 
which are not present in i?(Pi,P 2 ), 



N 



T c (Pi,P 2 )=^s i (P 1 )s i (P 2 ). 



(9) 



To obtain a normalized 2-particle correlation C(Pi,P 2 ), 
we choose the normalization 



iV(Pi,P 2 ) = KPi) i/(p 2 ) - T c (Pi, P 2 ) . 

i?(pl,p 2 ) 



C(Px,P 2 ) 



iV(Pi,P 2 ) 
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This choice is motivated by the experimental praxis 
of "normalization by mixed pairs" : An uncorrelated 
(mixed) pair is described by an unsymmetrized product 
state 



(X,Y,t)=/|-)(X,t)/(;)(Y,t), (11) 



for which the two particle Wigner phase space 
density and the corresponding detection probability 
PH ncorr (Pi,P 2 ) ca n be calculated according to (§). Tak- 
ing both distinguishable states $^ ncorr and $|? corr into 
account, and summing over all pairs we find 



A(Pi,P 2 ) = ^7 corr (Pi,P 2 ). 



(12) 



Hence, the normalization (10a) is the 2-particle detection 
probability for uncorrelated pairs. The correlator reads 
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C(Px,P 2 ) = l 



EL^( p i 5 p 2)e 
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i/(Pi) i/(P 2 ) - T c 



(13) 



In contrast to (g), this is a continuous function of the 
measured momenta Pi, P 2 , i.e., no binning of the cor- 
relator is required. Note that the normalization (10a) 
ensures that the correlator (|l3|) is always smaller than 2 
and equals 2 for Pi — P2 = 0. (This follows from the 
fact that the sum of the first two terms in (|7a| ) is always 
larger than the third one.) For any boson source, defined 
by an arbitrary phase space distribution S and a spa- 
tial wavepacket width er, Eq. ( [i"3"| ) provides an algorithm 
of how to calculate the 2-particle correlator, using (ffb|), 
©, © and (§). 

To understand how the correlator (|l3|) relates to the 
result of the coherent state formalism (1), the limit of a 
large number N of emission points is relevant. T c (Pi, P 2 ) 
in ( p"3| ) is a sum over N terms while the other terms in 
the nominator and denominator are sums of N 2 terms. 
In this sense, the Independence of (|l^) provides a finite 
multiplicity correction and can be neglected as a sublead- 
ing -i-contribution in the large iV limit of (jT" 
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lim C(Pi,P 2 ) = 1 

TV— >oo 



E£i«i(K)e<*' 



(E,=i«i(Pi)) E^i*i(Pa) 



(14) 

We note that in the derivation of (1), subleading -h- 
contributions are dropped Qj. The large N approxima- 
tions (1) and (|l4|) are clearly justified for pion interfer- 
ometry in ultrarelativistic (Pb-Pb) heavy ion collisions 
where typical pion multiplicities are in the hundreds. For 
smaller systems, however, and especially in studies of the 
multiplicity dependence of HBT correlations [15|, one 
might wish to start from the expression for finite mul- 
tiplicity (|l3"|). Expression (^) can be obtained from the 
coherent state result (1) by inserting 



S(x,K)=J2 s i( x > K )> ( 15a ) 

i=l 

Si(x, K) =J\fS(t — U) e~~^( x _ 

xe -a 2 (K- Pi ) 2 5 (15b) 



where Af is an arbitrary normalization factor. In this 
sense, (15a) is the emission function for a source £ with 
initial spatial localization a. It contains the informa- 
tion about how the initial phase space emission points 
Zi and the measured momenta K are correlated. Spa- 
tial and temporal components are not treated equally in 
( [15b] ), since our derivation is not Lorentz covariant. The 
Lorentz covariant setting used in (1) allows for an addi- 
tional dependence of S(x, K) on the temporal component 
of K which does not exist in our derivation. In practical 



applications however, the emission function (1) is used 
in the so-called on-shell approximation, where this addi- 
tional -FTo-dependence is not employed, ||. 

Both the 2-particle correlator (O) and the 1-particle 



spectrum v(P) in (8b) depend on the initial spatial local- 
ization a which is an additional free parameter. We now 
discuss this er-dependence. We first consider the limit 
(7 — > 0, in which the Gaussian wavepacket (|4a|) describes 
at freeze out (t — t"') a state with position uncertainty 
Ax = 0, i.e., the source is sharply ("classically") localized 
in configuration space. The prize for this optimal spatial 
information is that nothing can be said about the initial 
momenta pi at emission, the one-particle spectrum v{P) 
is flat. The measured momentum correlations contain 
spatial information about the source, namely 



lim C(Pi,P 2 ) 

a— >0 



1 + 



E (i ,,) cos((< Zi -g J ).(Pi-P 2 )) 



N(N- 1) 



(16) 



Due to the cos-term, the dependence of the 2-particle 
correlator ([l6]) on the measured relative momentum Pi — 
P 2 gives information on the initial relative distances qi — 
qj in the source. This is the HBT effect. Eq. (16) differs 
significantly from the cos-prescription (|J): here, Pi — P 2 
is the measured relative pair momentum, while p$ — pj 
in (||) denotes the initial momentum difference. As a 
consequence, the sum E(ij) m (0) S oes over a ^ pahs 
irrespective of the momenta pi, Pj since in the limit a — > 
0, all information about these initial momenta is lost, 
while the sum in (|^) goes only over those pairs for which 
the initial relative pair momentum pi —pj lies in the same 
bin as the measured Pi — P 2 . Since the correlator ( |l6| ) 
is a limiting case of (|i"4|), it is always larger than unity. 
In contrast, due to the wrong pair selection criterion, the 
correlator (|2|) can drop below 1. This insufficiency of 
(^|) becomes more significant for sources with strong q-p 
position- momentum correlation, as was noticed in Q. 
The other limiting case of (]l^) is the plane wave limit 



lim C(Pi,P 2 ) = l + <5 Pl! p 2 



(17) 



In this limit, nothing can be said about the spatio- 
temporal extension of the source since the 2-particle sym- 
metrized wave functions (|^) contain no space-time infor- 
mation. 

The difference between (0) and (|l7| ) shows that the 
er-dependence of the 2-particle correlator cannot be ne- 
glected. As pointed out already in | p)|JTTl , none of the 
two limits is realistic. For a — > 0, one has sharp infor- 
mation in configuration space but the momentum space 
information is lost and hence, the set S of phase space 
emission points (qi, Pi, tj) contains no information about 
the one-particle momentum spectrum is(P). In the limit 
a — > 00, on the other hand, no space-time information is 
contained in E. A realistic width a hence lies in between 
these two extremes. 
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For further discussing the a-dependence of 1- and 2- 
particle spectra, we now generalize Eq. ( |l5| ) to continu- 
ous phase space distributions p(q, p, t) which can encode 
statistical assumptions. To this aim, we extend the sum 



(15a) to a phase space integral weighted by a (classical) 
distribution p, [p|,|l3| 



S(x,K) 



J <i 3 q 4 d 3 pi dU p(m, pi, U) Si(x, K) . 



(18) 



The one-particle spectrum is obtained from ( |l8| ) via 
J d 4 x S(x, K). Especially, for a Boltzmann distribution 

2 

of temperature T, p(q,p, t) oc exp(— ^j), this one- 
particle spectrum has an effective temperature pT| 



T 



2mu 2 



(19) 



Hence, the one-particle spectrum broadens significantly 
for a narrow spatial width a. a is a free parameter which 
has to be determined from a comparison to data. How 
can this be done? One idea is to look at systems which 
can be expected to provide very small, almost pointlikc 
boson emission regions. Candidates are e.g. the p-p an- 
nihilation process jlij or Zo-decays [Q. The width of 
the HBT-correlator determined for these systems should 
be dominated by the width a. To obtain an argument 
supporting this idea, we consider the extreme case of a 
"pointlike source" p with no momentum dependence, for 
which all particles are emitted from the same space-time 
position q, t. Calculating the emission function (18) for 
the corresponding /o(qj, Pi, = 5( 3 '(qj — q)S(t — ti), we 
find 



C(Pi,P 2 ) = l + e-ff 



^Q 2 



(20a) 
(20b) 



Several assumptions enter this result: for pointlikc 
sources with an additional momentum dependence, the 
correlation is in general more complicated. Also, the 
width a could in principle depend on the emission points 
Zi, the localized wavepackets could have a different, non- 
Gaussian shape, etc. Still, Eq. ( po| ) suggests that the size 
of the HBT radius parameters measured for very small 
boson emitting systems is essentially given by a. 

From the pion interferometric measurements of sys- 
tems like the p-p annihilation process or Zp-decays 
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one infers on the basis of ( pOq ) a pion wavepackct 
width of the order a ~ 1 fm. This is in good agreement 
with the natural localization scale of the pion, its Comp- 
ton wavelength H] . Remarkably for such a localization, 
the additional quantum contribution to the temperature 
T off in @ is of order w 100 MeV. This indicates 

that the initial spatial localization width a plays an im- 
portant role in accounting for the slope of the measured 
one-particle spectra. 

In the present formalism, the role of an event generator 
for the boson emitting source is to provide a dynamical 



calculation of the phase space occupation E from some 
more fundamental initial condition. The current praxis 
for event generators of heavy ion collisions amounts to 
determining the 1-particle spectrum in the limit a — > oo. 
We have shown that this limit is unrealistic and that a 
realistic spatial wavepacket width leads to a substantial 
broadening of the one-particle spectrum. Our main re- 
sult is an algorithm which allows for the calculation of 
both the one-particle spectra via v(P) in (Be), and the 
two-particle correlations via C(Pi,P2) in (IS), starting 
from an arbitrary initial phase space distribution E of 
wavepackets with arbitrary spatial localization a. Here, 
the spatial width a is an additional free parameter which 
has to be fixed in comparison with experimental data. 
The slope of the transverse mass spectra (19) provide a 
lower bound for a while the 2-particle correlators pro- 
vide an upper bound. For pions these bounds are very 
tight and a realistic width a is of the order of the pion 
Compton wavelength. 
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